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We show that a simple rate-and-state theory accounts for most features of both time- independent 
and time-dependent plasticity in a spatially inhomogeneous situation, specifically a circular hole in 
a large stressed plate. Those features include linear viscoelastic flow at small applied stresses, strain 
hardening at larger stresses, and a dynamic transition to viscoplasticity at a yield stress. In the 
static limit, this theory predicts the existence of a plastic zone near the hole for some but not all 
■^- ' ranges of parameters. The rate-and-state theory also predicts dynamic failure modes that we believe 

may be relevant to fracture mechanics. 
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I. INTRODUCTION 



Since the work of Hart H in the 1960's, scientists have understood that a satisfactory theory of plastic deformation 
in solids must include dynamic variables that describe the internal states of materials. The deformation field itself 
cannot be sufficient. It cannot, in any natural way, describe the irreversible changes that lead to hysteretic stress- 
strain curves, or to the transition from nonlinear viscoelastic to viscoplastic behavior with increasing applied stress. 
Conventional theories of plasticity cope with these limitations by specifying phcnomcnological rules to suit various 
situations and histories of deformation. For example, strain hardening curves, viscoplastic laws, or the distinctions 
between loading and unloading behaviors are determined from experiment and used as needed in computations. In 
most treatments there is also a sharp distinction between time-independent and time-dependent behaviors, with little 
O or no indication of how these properties may be related to one another. We believe that a deeper, more nearly 
fundamental level of phenomenology is required for modern applications, for example, for computing deformations 
,_^ ' near moving crack tips. Plasticity is an intrinsically time dependent phenomenon; time-independent descriptions 
£> . should emerge as static limits of fully dynamic theories. 

In a recent paper J2|, M. Falk and one of the present authors (JSL) proposed a theory of plastic deformation in 
amorphous solids in which they introduced an internal state variable to describe the orientations of what they called 
"shear-transformation zones." (We refer to this as the "STZ" theory.) The resulting "rate-and-state" equations (a 
concept that is widely used in the seismological literature [SEf and in recent theories of friction ||) successfully 
describe the full range of viscoelastic and viscoplastic phenomena, including hysteretic effects. The basic structure 
of the STZ theory appears to be broadly applicable. The new internal state variable might equally well describe, for 
1*>! . example, anisotropy in the way dislocations pile up near defects in crystalline materials. 

Reference [|2| discusses plasticity only in spatially uniform situations. Our purpose here is to apply a simple version 
of the STZ theory to a spatially inhomogeneous situation and to make contact with conventional plasticity theory. 
We especially want to learn whether the conventional picture of a time-independent plastic zone appears in the static 
limit. Looking ahead to fully dynamic situations such as fracture, we also want to understand the dynamics of plastic 
flow in regions of concentrated stress. We shall show that the conventional time-independent concepts — the "yield- 
surface" hypotheses — do emerge from dynamic theories in an approximate way in many normal situations. As we 
shall argue, however, the rate-and-state theory is simpler, richer and more general than the conventional approaches. 

The problem of describing spatially inhomogeneous plastic deformation is best approached by looking at a simple 
example where questions of time dependence and compatibility are not obscured by mathematical details. A growing 
circular hole in a very large stressed plate satisfies the criterion of simplicity. It shares important features with the case 
of plastic deformation near a crack tip; the tractions are applied at a distance and the stresses are concentrated near 
the hole. Several researchers have addressed the problem of dynamic hole growth in the context of ductile fracture H 
and spallation M . Carroll and Holt and later Johnson included inertial effects but neglected rate dependent plasticity 
p,M . Bodner and Partom used a version of a conventional rate dependent plasticity theory but did not study a stress 
controlled situation M. 

The scheme of this paper is as follows. In Section O, we derive equations of motion for the radius of the circular hole 
and the surrounding stress field, assuming a general form for the constitutive relation that governs the rate of plastic 



deformation. Section III contains a brief summary of several results of conventional plasticity theory which will serve 



as points of comparison for the rate-and-state analysis. A simplified version of the STZ model is introduced in Section 



IV A where, for completeness, we outline some important properties of this model that were repor ted p reviously in 
]2| and []10|. We describe the dynamic behavior of the STZ model for the hole problem in Section IV B . The paper 
concludes with remarks about the implications of these results in Section M. 



II. DYNAMIC PLASTICITY IN A CIRCULAR GEOMETRY 



Throughout this analysis, we consider only a two-dimensional solid in a state of plane strain, and assume that 
inertial effects are negligible. Suppose that a circular hole in this system has radius R(t) at time t. Outward tractions 
at the distant edges of the plate cause the pressure p far from the hole to be p — * — a^. We introduce polar coordinates 
r and 9 that define an Eulerian reference frame, so that u(r, t) is the radial displacement, measured from some initial 
reference state, of the material currently at position r. The function u(r,t) is the only degree of freedom in the 
problem. 

The total rate-of-deformation tensor (including both elastic and plastic parts) is diagonal with components: 



T) tot = —■ 

rr g ■ 



TV 



where v is the material velocity: 



du/dt 
1 — du/dr 



(2.1) 



(2.2) 



For small strains, the tensor V is approximately equal to the total strain-rate tensor e . Eq.(2.1) implies that the 
components of V tot satisfy the compatibility condition: 







-(rVl°i)=Vl° r K 



(2.3) 



The stress tensor has the form: 

drr = -P - s, aee =—p+ s. (2.4) 

Here, s = (agg — a rr )j1 is the deviatoric stress. The elastic stress-strain relations are: 

2/*e£ = -(l-2i/)p-a, 2(j,sf e = -(1 - 2u)p+ a, (2.5) 

where /i is the shear modulus and v is Poisson's ratio. In the absence of inertial effects, balance of forces implies 



dp 1 d 2 

or r z or 

The boundary condition at r = R is 

- <7 rr (R, t) = p(R, t) + s(R, t) = 0. 
For simplicity, we neglect surface tension. Also at r = R, we have 

R 



R 



= -Dlf[ S (R),...}. 



Here and elsewhere, dots above symbols denote derivatives with respect to time t. 
For incompressible plasticity, the constitutive relation has the form: 



(2.6) 



(2.7) 



(2i 



(2.9) 



where T> is some function of s and possibly other variables as indicated by the ellipsis. Combining (2.9) with the 
compatibility condition (2.2), the elasticity equations (2.5), and force balance (|2.6D, we find 



©(*,...) + - -S= -At 1 . 
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(2.10) 



wher e C(t) is an r-independent constant of integration. Further analysis using the boundary conditions (2.7) and 
(U) at r = R yields 



C(t) = 



and 



l-2u 



s(R,t) RR; 



(2.11) 
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(2.12) 



Eqs. ( |2. ld|) through ( 2.12 ), supplemented by equations of motion for other arguments oiV(s, ...), constitute a coupled 
set of first-order differential equations suitable for computing the time evolution of s(r,t) and R(t). 

If R(t) is the only length scale in the problem, we can look for self-similar solutions in which R/R = uj = constant 
so that the hole radius i s gro wing exponentially. All functions depend only on £ — r/R{t). Combining the expression 
for C{t) in ( 2.11 ) with ( |2.10 ), and transforming to the scaling variable £, we find 
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Force balance (2.6) plus the boundary condition (2.7) imply 



If these self-similar solutions exist for u> > 0, they describe unbounded plastic failure of the material. 



(2.13) 



(2.14) 



III. CONVENTIONAL THEORIES 



In a typical time- independent approach to this problem JETl, one assumes that there exists a maximum value of 
s, say s y , and that, in any sufficiently slow deformation, the material adjusts its state so that s < s y everywhere. 
Technically, the condition s = s y is a special case of a Tresca yield surface in the space of stress components. For 
present purposes, we take this assumpt ion to mean that the hole is surrounded by a plastic zone, R < r < R\, 
within which the force-balance equation ( |2.6| ) remains valid but the condition s = s y replaces Hookean elasticity (|2.ED . 

Outside this zone, r > R%, p = — a^ and s — s y R\/r 2 . Continuity of stress at R\ means that, within the zone, 

P 



(Too — 2s y ln(r/i?i). Then the boundary condition (2.7) at r — R implies that 
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(3.1) 



Thus, these assumptions predict that a stationary state with a non- vanishing plastic zone exists for a^ > s y . A 
calculation of the displacements similar to that described by Hill |]12) shows that the hole radius R diverges as a^ 
approaches an upper threshold stress c^ which (for the case s y j\x <C 1) is given by: 
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(3.2) 



To see what happens above this stress, we must consider a time-dependent theory 

The simplest conventional time-dependent hypothesis is the Bingham law which, for s > 0, we can write in the 
form: 



V{s) = 



_ J a(s — s y ) for s > 







otherwise, 



(3.3) 



where a is a response coefficient. Essentially by definition, the combination of (^3) with ( 2.10J )— ( 2. 12 ) describes an 
elastic perfectly-plastic material whose stationary states, for s y < a^ < tr^, are the sam e as those des cribe d in th e 
preceding paragraph. We have checked that these states are stable attractors by using (3.3) to integrate ( 2.10| )-( f2.12 ), 
usin g the initial condition s(r, 0) = a oo R 2 (0)/r 2 (for r > R(Q)). We found no surprises; a plastic zone consistent with 
(EO) forms around the growing hole. 



We also can compute the self-similar solutions of (2.13) for the Bingham model. For ease of analysis, we write these 
for the case of incompressible elasticity, v = 1/2. In the outer, elastic region, £ > £1 = (fi/sy) 1 ^ 2 : 



m 



a 



and, in the plastic region, s > S v , 1 < £ < £i: 



m 
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jU,a + cj / £ 
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(3.4) 



(3.5) 



where f3(w) = 2/j,a/u>. Note that the exponent /3 becomes indefinitely large in the limit of small u>; thus the second 
term in the square brackets in ( pJ5|) produces a function s(£) that is sharp ly bent but continuous at £ = £i. To 
complete the calculation, we use ( |3.5D to evaluate the right-hand side of fl2.14| ). After rearranging and taking the limit 
of small w, we find: 






(3.6) 



where tr^ is the upper threshold defined in (3.2). Thus, as expected, the dynamic failure modes start where the 
time-independent theory breaks down. 



IV. THE STZ MODEL 



A. Basic properties 



For present purposes, it will be sufficient to use the simplified version of the STZ model that we introduced in an 
ear lier one-dimensional analysis ftO] . As in the Bingham case, the model is specified by the function T>(s, ...) defined 
in (2.9). For two-dimensions, with circular symmetry, we write: 



V(s,A) = -(Xs-A) 



(4.1) 



and supplement this by an equation of motion for the state variable A: 

A = X>(s,A)(l-7sA). 



(4.2) 



A(r, t) is the single independent element of a diagonal, traceless tensor which describes the local anisotropy of the 
shear-transformation zones. We are omitting the other state variable in H that describes the density of STZ's on 
the assumption that this quantity quickly reaches its equilibrium value. More importantly, this simplified version of 
the STZ model omits the strongly s-dependent rate factor that governs memory effects. This version is qualitatively 
sensible if we load the system only once in one direction, but it does not behave properly if the loading is cycled in 
any way. 

The inverse stress 7 can be eliminated in favor of a group of parameters that plays the role of a dynamic yield 
stress, specifically, s y = l/\f\rj. Note what is happening here for spatially uniform situations. For s < s y , A(i) has 
its stable fixed points on the line A = As, where T> = e p ge — 0. In this region, the material is nonlinearly viscoelastic. 
For s <C s y or, equivalently, 7 — > 0, it obeys a conventional creep-compliance law flO(| : 



t (t) = (l + X)s{t)-X 



/dt' exp 
-00 



1 



(t-f) 



s(t'). 



(4.3) 



We obtain a particularly important result by supposing that the system is initially in a state with e v ss = A = and 
that a stress s < s y is applied suddenly at time t = 0. A simple calculation then yields: 
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(4.4) 
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e p (t) approaches £g na j exponentially in time with a relaxation time T rc i ax that diverges as s — ► s y 
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1 - {s/Syf 



(4.5) 



Eq. (4.4) is a strain-hardening curve; that is, £^ nal is the non-recoverable plastic strain produced after an infinitely 



long time by the deviatoric stress s. In the limit A — > 0, with s y held constant, 



e finai vanishes for s 



< s, 



but 



have any positive value for s = s y . Thus the parameter A is a measure of the deviation from perfect plasticity. The 

diverging relaxation time near s — s y , however, has no simple analog in conventional descriptions of strain hardening. 

This one-parameter fit (4.4) to the shape of the strain-hardening curve is much too simple even for the fully 



nonlinear STZ model, where the behaviors at small stresses and at stresses near s y are determined by different groups 
of parameters. Moreover, the small-s behavior of this t — > oo curve is not what is measured experimentally. In real 
materials and in the full STZ theory, the plastic deformation rates at small stresses are too small to be observed, and 
the material behaves as if it were purely elastic. For an illustration of this behavior, see Fig. 5 in J2J]. 
For s > s y , A goes to 1/js and 



TS 



)=—(«- Sy). 



V 



(4.6) 



In dynamic situations at large stress, therefore, the STZ model looks like a Bingham plastic. In short, even this highly 
simplified version of the STZ model describes much of both static and dynamic plasticity. 



B. Dynamically growing hole 



We return n ow to the circle problem. As a first investigation, we have numerically integrated (2.10)— (2.12), sup- 
plemented by (4.2), to find the time evolution of s(r,t), A(r,t), and R(t). In Figs. [I] and ||, we show what happens 



if we suddenly apply the stress s(r, 0) = a c 
previously undeformed system, A(r, 0) = 0. 



R 2 (0)/r 



We set Xs y = 0.005. 



Sy = 0.1/i, and <7oo = 0.2/x and assume a 



For these values of the parameters, the hole grows for a while and then 



stops, and a plastic zone with s = s y , A = Xs y forms around it. Apart from the fact that the outer boundary of the 
plastic zone is smooth rather than sharply defined, this static limit of the time-dependent deformation is qualitatively 
consistent with conventional, time-independent plasticity theory. 

In Fig. ||, we show an analogous set of curves for a substantially larger value of A, specifically, Xs y — 0.5. According 
to (4.4), this system deviates appreciably from perfect plasticity. A plastic zone does form around the hole, but it 
has no sharp outer boundary. If we estimate the position of this boundary, say, by finding the point of inflection in 
the final curve s(r), we find that the relation (3.1) is strongly violated. 

As in the Bingham model discussed in Section III , the STZ model predicts the existence of unbounded failure modes 
for sufficiently large a^. The stationary states of the kind shown in Figs. H and @ cease to exist beyond a^. Indeed, 
as we show in Fi g. El) the ra dius of the hole diverges at a^. We find these dynamically growing states via the scaling 
analysis of Eqs. (2.13) and ( |2.14| ). It is useful to make the following changes of variables: 



R 2 (t) 



- = C; s(r, t) = 8y il>{C)\ A(r, t) = As^(C). 



(4.7) 



We find: 



2(1 — u)ujt _ dib lot _ 

^-^ + - r^ — C-7^ = ^ — C 

/iA ac, \s y 



1 + (1 - 2v) 2^(1) 



The equation of motion for A, i.e. (4.2), becomes 



2^tC— = (ip-<p)(l- <pip). 



Finally, ( 2.14 ) becomes: 



c 



HO- 



(4.8) 



(4.9) 



(4.10) 



As a first step in interpreting these equations, we compute the threshold a^ by taking the limit u> — > 0. To avoid 
unnecessary complication, we again set v = 1/2. In this case, ( [2.13 ) and (4.S) reduce to f/ 1 ~ f and, after a simple 
integration: 



Xsy In 



i + y> 

l-t/j 






Note that, if s y /fi<^l 1 (which is generally true for realistic situations), then 



ip(() = tanh 



( 



2As, 



(4.11) 



(4.12) 



This solution exhibits a plastic zone with a smooth elastic-plastic boundary only for Xs y <C 1, in which case there is 
a region between £ = 2\s y and £ = 1 in which ^»1, 

It is easy to compute "0(C) without making the latter approximation and, via (4.1CJ), to obtain the threshold stress 



■rth 



If (s y /n)(l + 2\[i) < 1, then 



£~ wi + lnfii) _ hu i -l 2A/M. 
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(4.13) 



which agrees with (3.2) when A = and v = 1/2, i.e. in the limit of perfect plasticity. Note that we have chosen the 
parameters in Fig. Ill to lie within this range. If, on the other hand, \s y ^> 1, then a^ rs 1/2A. Here the threshold lies 
below s y ; the material is highly deformable and conventional plasticity theory has no range of validity. We illustrate 
this point in Fig. a. 

Expanding the solutions of (4.8) and (4.9) to first order in w, 

Sy/ ' \l <C 2XSy <C 1, 



we find the following behavior near threshold: For 



2A 



[1 + 2Xs y ln(As y )] (<7 
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(4.14) 



and, for Xs y ^$> 1, 



- (ct c 
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(4.15) 



In each of the last two equations, the quantity a^ 
paragraph. 



has the value computed in the corresponding limit in the previous 



We have studied the behavior of the similarity solutions (4.8)-(4.1C) numerically. Just as in the Bingham plastic, 



the deviatoric stress at the surface of the hole in the STZ material grows with the hole expansion rate u as shown 
in Fig. @. While intuitively obvious, this result is relevant to understanding stress transmission to brittle crack tips. 
In Fig. we show the A-dependence of the stress for a slowly expanding hole. This solution is essentially the u> — ► 
limit obtained in Eq. (4.11). In a "softer" material with larger A, the plastic zone shrinks and disappears completely 
for a large enough A. This softening of the material for larger A leads to the decrease in the threshold stress c^ above 
which dynamic failure modes exist. 



V. DISCUSSION 



Our principal point is that the STZ model, a simple example of a rate-and-state theory, provides an extremely 
compact and physically motivated description of essentially all of plasticity theory both static and time dependent. 



In just two constitutive relations, (4.1) and (4.2), containing just two dimensionless groups of parameters, s v j \x and 



Xs y , plus one time constant r, we capture linear viscoelasticity at small stress, strain hardening at larger stress, and 
a dynamic transition to viscoplasticity at a yield stress. All of these properties have been described previously in |2| 
for homogeneous situations and in [lQ ] for an inhomogeneous one dimensional case. 

In this paper, our principal interest has been to make contact with conventional theories of plasticity by looking 
at deformation near a circular hole. As a rule, we recover conventional results when both s y <C fi and Xs y -C 1. The 
first of these conditions, that the yield stress must be much less than the shear modulus, is generally true for ordinary 
materials. The second condition, according to (4.4), says that the plastic strain at which large-scale deformations 
start to occur must be much smaller than unity. This too is ordinarily satisfied by rigid solids. However, it remains to 
be seen whether this A-dependent condition might be modified in more realistic parameterizations of rate-and-state 
theories. 

One place where the rate-and -sta te theory differs qualitatively from conventional plasticity is in the interpretation 
of the strain hardening curve (4.4). Here, this curve must be interpreted as the t — > 00 limit of the dynamic 



inelastic response to an externally applied stress, not as an instantaneous response that is mathematically equ ival ent 



to nonlinear elasticity. Moreover, the parameter A that dete rmines the deviation from perfect plasticity in (4.4) is 



the same parameter that appears in the viscoplastic law (4.6), and is also the same parameter that determines both 
the shape of the plastic zone and the threshold a 1 ^ at which static solutions give way to time-dependent, unbounded 
failure modes. 

Perhaps the single most important advantage of the STZ theory is that the material in the plastic zone is charac- 
terized not just by the stress and corresponding displacement fields but also by the state variable A. This variable 
tells us in a natural way that the material inside the zone will respond differently to subsequent changes in stress than 
will the undeformed material outside it. To take advantage of this feature, however, we must go beyond the truncated 
version of the STZ model that we have used here and include the strongly nonlinear, s-dependent rate factor derived 

in i- 

Without doing any further calculations, we can sec how the nonlinear rate factor produces memory effects. Suppose 
that we start in a stressed configuration such as the one shown in Figs. [I] and in which a plastic zone has formed 
around the hole; and suppose that we then unload the system by quickly reducing the external stress to zero. At s = 0, 
the rate factor becomes extremely small; in the absence of a stress that can induce transitions from one orientation 
to another, the zones remain as they were in the previously stressed state. That is, A(r) remains unchanged on 
unloading. The system relaxes elastically, but the plastic degrees of freedom do not revert to their earlier values. 
Accordingly, there must be residual stresses in the region near the hole. The material will "remember" the magnitude 
and direction of its previous loading because it "knows" the function A(r), which will determine via the nonlinear 



generalizations of the constitutive equations (4.1) and (4.2) what happens when new stresses are applied 



The other feature that we have emphasized in this analysis is the dynamic failure that occurs at aj^ , a feature that 
also occurs in the conventional time-dependent theories. We believe that the dynamic growth modes may provide a 
clue for solving a long-standing puzzle in the theory of fracture, specifically, the question of how breaking stresses can 
be transmitted through plastic zones to crack tips. This puzzle pertains to quasi-static crack advance, which would 
seem to be impossible if, as in conventional plasticity theory, the stresses near a crack tip are constrained to be less 
than or equal to the yield stress. The new modes, especially those that occur when there is appreciable deviation 
from perfect plasticity, raise the possibility that, like the growing hole, crack growth could occur via plastic flow near 
the tip. We hope to explore this possibility in a subsequent report. 
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FIG. 1. Deviatoric stress s/ fi as a function of the distance from the center of the hole shown at several times after the 
application of the stress (Too = 0.2p. The yield stress s y — 0.1/x, Poisson's ratio v = 0.3, Xs y = 0.005. 
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FIG. 2. Time evolution of A for the same values of the parameters as in Fig. 
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FIG. 3. Same as in Fig. hi but wuth Xs y = 0.05. Note that the approach to a stationary state is faster and that the plastic 
zone is much less pronounced. 
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FIG. 4. The final radius of the hole versus the applied stress in units of the corresponding threshold stress for three values 
of Xs y . s y — 0.1/i. The solid line is the prediction for the Tresca plastic. 
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FIG. 5. Threshold stress in units of the yield stress as a function of A. Note that for soft materials, the thresold for unbounded 
failure is smaller than the yield stress. 
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FIG. 6. Similarity solutions s(r/R(t)) for three different values of the hole growth rate u>. \s y — 0.005, s y = 0.1/i. 
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FIG. 7. Similarity solutions s(r/R(t)) for three different values of Xs y in the case of a small growth rate lot = 0.0001. 
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